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We propose a new method to determine magnetic fields, by using the magnetic-field in¬ 
duced electric dipole transition ^P 3/2 in Fe®'*' ions. This ion has a high 

abundance in astrophysical plasma and is therefore well-suited for direct measurements of 
even rather weak fields in e.g. solar fiares. This transition is induced by an external magnetic 
field and its rate is proportional to the square of the magnetic field strength. We present 
theoretical values for what we will label the reduced rate and propose that the critical energy 
difference between the upper level in this transition and the close to degenerate 3p^ id ^ 05/2 
should be measured experimentally since it is required to determine the relative intensity of 
this magnetic line for different magnetic fields. 


I. INTRODUCTION 

One of the underlying causes behind solar events, such as solar flares, is the conversion of 
magnetic to thermal energy. It is therefore vital to be able to measure the magnetic field of the 
corona, over hot active areas of the sun which exhibits relatively strong magnetic fields. In order to 
follow the evolution of a solar flare, continuous observations are required, either from space or by 
using a network of ground-based instruments. It is therefore unfortunate that there are no space- 
based coronal magnetic field measurements, but only model estimates based on extrapolations 
from measurements of the photospheric fields (Schrijver et al. 2008). Ground-based measurements 
are performed either in the radio range (White 2004) across the solar corona, or in the infrared 
wavelength range (Lin et al. 2004) on the solar limb. Infrared measurements of magnetic fields 
are limited by the fact that the spectral lines under investigation are optically thin. On the other 
hand, gyroresonance emission is optically thick, but refers only to a specific portion of the corona, 
which has a depth of around 100 km. From these measurements an absolute field strength at the 
base of the corona, above active regions, in the range of 0.02 - 0.2 T was obtained (White 1997). 

In this work we present a completely new method to measure magnetic fields of the active corona. 
This method is based on an exotic category of light generation, fed by the plasma magnetic field, 
external to the ions, in contrast to the internal fields generated by the bound electrons. The 
procedure relies on radiation in the soft x-ray region of the spectrum, implying a space based 
method. This magnetic-field induced radiation originates from atomic transitions where the lifetime 
of the upper energy level is sensitive to the local, external magnetic field (Beiersdorfer et al. 2003, 
Li et al. 2013, Grumer et al. 2013, Grumer et al. 2014, Li et al. 2014). We will show that there is 
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Figure 1: (Color online) Schematic energy-level diagram for Chlorine-like ions with Z < 26 and zero nuclear 
spin, where ^ 07/2 is the lowest level in the configuration 3s^3p^3d. For ions with Z > 26, a level crossing has 
occured and ^ 05/2 is lower than ^ 07 / 2 - Under the influence of an external magnetic field, an El transition 
opens up from the ^ 07/2 to the ground state through mixing with the ^ 05 / 2 - 


a unique case where even relatively small external magnetic fields can have a striking effect on the 
ion, leading to resonant magnetic-field induced light, due to what is called accidental degeneracy 
of quantum states. 

The impact of the coronal magnetic field on the ion is usually very small due to the relative 
weakness of these fields in comparison to the strong internal fields of the ions. The effect therefore 
usually only contributes very weak lines that are impossible to observe. However, sometimes the 
quantum states end up very close to each other in energy, they are accidentally degenerate, and 
the perturbation by the external field will be enhanced. If this occurs with a state that without 
the field has no, or only very weak, electromagnetic transitions to a lower state, a new and distinct 
feature in the spectrum from the ion will appear — a new strong line. Unfortunately, since the 
magnetic fields internal to the ion and externally generated in the coronal plasma differs by about 
five to seven orders of magnitude, the probability of a close-enough degeneracy is small. But in 
this report we will discuss a striking case of accidental degeneracy in an important ion for studies 
of the sun and other stars, Fe®"*". 

The origin of the new lines in the spectra of ions is the breaking of the atomic symmetry by 
the external field, which will mix atomic states with the same magnetic quantum number and 
parity. This will in turn introduce new decay channels from excited states (Andrew et al. 1967, 
Wood et al. 1968), which we will label magnetic-field induced transitions (MITs) (Grumer et al. 
2014). These transitions have attracted attention recently, when accurate and systematic methods 
to calculate their rates have been developed (Grumer et al. 2013, Li et al. 2013). 
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II. STRUCTURE OF CHLORINE-LIKE IONS AND MITS 

The structure of the lowest levels of Chlorine-like ions is illustrated in Figure [H The important 
levels in the present study are the two lowest in the term 3p^3d which turn out to have very 
different decay modes. Withont external helds and ignoring the effects of the nuclear spin, they 
both decay to the 3p^ ^ f^ie gronnd configuration, but while the J = bl2 has a fast 

electric dipole (El) channel, the J = 112 can only decay via a slow magnetic qnadrupole (M2) 
transition. In the presence of an external magnetic field, these two states will mix and induce an 
El, MIT competing transition channel from the J = 7/2 level. Eor most ions the M2 transition 
is still the dominant decay channel, bnt a crossing of the fine structnre levels '^ 07/2 and ^ 05/2 
between Cobalt and Iron (see Eigure. [2|) will change the pictnre. As a matter of fact, for Iron this 
fine-structure splitting energy is predicted to be at a minimum and the MIT-contribution to the 
decay of the J = 7/2 level will be strongly enhanced. 



Element (Z) 


Figure 2: The fine-structure splitting energy AE = ^(^ 07 / 2 ) — as a function of the nuclear charge, 

along the isoelectronic sequence from calculations reported here. The dashed line in green marks AE = 0. 

Unfortunately, there are large variations in the predicted value of this energy difference (see 
Table (H) . The aim of this report is therefore to (a) make a carefnl theoretical study of the energy 
splitting between the two '^D'^-levels along the isoelectronic sequence, to confirm the close degener¬ 
acy for iron, and (b) to make an accurate prediction of what we will label the reduced decay rate, 
next section). This reduced rate can be combined by the experimentally determined 
wavelength and energy splitting to give the MIT-rate for different magnetic fields. 
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Table 1: ATOMIC DATA FOR FE X 


Method 

A 

AE 

Aei 

Observation Solar (Thomas et al. 1994, Brosius et al. 1998) 

257.25 

0 


Solar (Sandlin 1979) 


5 


present 

257.7285 

20.14 

6.30[6] 

MCDF (Huang et al. 1983) 

246.4924 

78 

1.63 [6] 

MCDF (Dong et al. 1999) 

256.674 

108 

6.27[6] 

Theory MCDF (Aggarwal et al. 2004) 


54.85 


MR-MBPT (Yasuyuki et al. 2010) 

257.1924 

18 


R-matrix (Del Zanna et al. 2012) 

246.8890 

109.74 


Cl (Bhatia et al. 1995) 

256.1974 

-58 

1.21[6] 

Cl (Deb et al. 2002) 

257.0846 

21 

2.42 [5] 


NOTES. -List of the wavelength A (in A) of the magnetically induced 3p'*3ci "^ 07/2 — t 3p® ^P 3/2 transition, the fine 
structure splitting AE = E ("^ 05 / 2 ) — E ('* 07 / 2 ) (in cm“^) and the rate A_Bi(3p'*3d '* 05/2 ^ 3p® ^P 3 / 2 ) of the El 
transition (in s“^) for Chlorine-like Iron. The second column lists the various sources. Numbers in square brackets 
are the power of 10. 


III. THEORETICAL METHOD AND COMPUTATIONAL MODEL 


A. Theoretical method 


The basis of our theoretical approach is described in our earlier papers on MITs (Li et ah 2013, 
Grumer et al. 2013). In our example the reference state is \3p^3d “^ 07 / 2 ), which we can represent 
by a mixture of two pure states in the presence of a magnetic field 

|“3/3d '^D 7 / 2 ” M) = do\3p'^3d ^ 07/2 M) + di{M)\3p^3d ^ 05/2 M). (1) 


where we ignore interactions with other atomic states, since their energies are far from the reference 
state. The total El transition rate from a specihc M sublevel in the mixed “3p^3d “^ 07 / 2 ” to all 
the M' sublevels of the ground level 3p^ ^P 3/2 can be expressed as: 


Amit{M) = 'Y2, M') 

M' 


2.02613 X 10^® 
3A3 


d,{M){3p^ 2p^^^||p(i)||3^43^ 


( 2 ) 


where A is the transition wavelength and di{M) depends on the magnetic quantum number M of 
the sublevels belonging to the 3p^3d ^ 07/2 level. For the 3p‘^3d '^ 07/2 level, di{M) is given by 


dx{M) 


(4D5/2M|ih^|4D7/2M) 
E(4D7/2) - E(4D5/2) 


-B 


49 - 4M2 (4D5/2||N(i) + AN(i)||4D7/2) 


63 


E(4D7/2) - E(4D5/2) 


( 3 ) 


As a result, the total rates of the 3p^3d ^^ 07/2 ^ 3p^ ^P 3/2 MITs from individual sublevels can be 
expressed as 


Amit{M) 


“mzt(M) 


A3(AE)2’ 


( 4 ) 


where B is in units of T, A is in units of A, AE = E(^Dy/ 2 ) ~ units of cm and we 

have defined a reduced transition rate as 


2.02613 X 10^® • (49 - 4M2) 


(4D5/2||NW + ANP)|| 4D7/2)(®P3/2||P^'^I|'D5/2) 




189 


( 5 ) 
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The reduced rate defined in this equation is independent of the transition wavelength, the magnetic 
field strength as well as the energy splitting. This gives us the property that relates the MIT-rates 
to the external magnetic field strength. To determine AmitIM), we recommend to use theoretical 
values of a^jrp{M), as reported in this work, combined with experimental values of the energy 
splitting and wavelength. 


B. Correlation Model 

The calculations are based on the Multiconfiguration Dirac-Hartree-Fock (MCDHF) method, in 
the form of the latest version of the GRASP2K program (Jonsson et al. 2013). A single reference 
configuration model is adopted for the even-parity(3p^3(i) and odd-parity(3p^) states, and the Is, 
2s, 2p core subshells are kept closed. The set of CSFs is obtained by single and double excitations 
from the n=3 shell of the reference configurations to the active set. The active set is augmented 
layer by layer to n=7 {Imax = 4) when satisfactory convergence is achieved. For each step, we 
optimize only the orbitals in the last added correlation layer at the time. In the final calculations, 
the total number of CSFs was 16490 for the odd-parity(J=3/2, 1/2) and 523421 for the even-parity 
(J = 1/2, 3/2, 5/2, 7/2, 9/2) cases. 

The resulting excitation energies of the ^ 07/2 and ^ 05/2 changes by less than 0.1% in the last 
step of the calculation. The final excitation energies agree with experiment to within 1%. The 
crucial energy splitting between “^ 07/2 and ^ 05/2 is well-converged, except for iron where the close 
degeneracy occurs. To better represent this critical case we extended the calculations to include 
single excitations from the 2s and 2p subshells. 


IV. RESULTS AND DISCUSSION 
A. Isoelectronic Sequence 

We present in Table [2] all the important properties, according to Eq. ([1|), involved in computing 
the MIT-rates, i.e. the reduced transition rate a^jrp[M), the energy splitting, AE, between the 
two levels 3p^3d ^05/2 and ^07/2, together with the wavelength. A, of the 3p^3d ^07/2 —> 3p^ ^P3/2 
transition. 

In the absence of an external magnetic field, magnetic quadrupole (M2) is the dominant decay 
channel for the 3p‘^3d ^07/2 — > 3p^ ^P3/2 transition. When an external magnetic field is intro¬ 
duced, an additional decay channel is opened and we define a average transition rate Amit of the 
3p^3d ^D 7/2 —>■ 3p^ ^P 3/2 transition. 


Amit 


Yjm 

2J+1 


( 6 ) 


Rates for any field can be obtained by eq. 0] and [5l We plot the transition rates 
A = Am 2 + Amit along the isoelectrionic sequence in Figure [3] (a) for some magnetic-field 
strengths and AE = 20.14 cm“^ for Iron. It is clear that the magnetic field influences the tran¬ 
sition rate substantially for Iron due to the close degeneracy. To further illustrate the resonance 
behaviour of this effect, we also used an astrophysical value (Sandlin 1979) AE = 5 cm~^ for 
Iron, in Figure [3] (b). 
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Table 2: CALCULATIONAL RESULTS FOR CL-LIKE IONS 


ions 





AE 

A 

Am2 

M = ±1/2 

M = ±3/2 

M = ±5/2 

Ar+ 

1.26[0] 

7.994[0] 

6.662[0] 

3.997[0] 

165.5 

762.5978 

K2+ 

3.20[0] 

7.686[0] 

6.405[0] 

3.843 [0] 

202.87 

600.4231 

Ca3+ 

6.16[0] 

7.621 [0] 

6.351[0] 

3.810[0] 

246.65 

500.1306 

Sc4+ 

1.03[1] 

7.910[0] 

6.592[0] 

3.955[0] 

283.9 

430.6129 

TiS-t 

1.59[1] 

8.443[0] 

7.036[0] 

4.222[0] 

307.02 

379.0310 

V6+ 

2.31[1] 

9.093[0] 

7.577[0] 

4.546 [0] 

306.35 

338.9583 

Cr^+ 

3.22[1] 

9.807[0] 

8.172[0] 

4.903[0] 

270.41 

306.7765 

Mn®+ 

4.33[1] 

1.055[1] 

8.793[0] 

5.276[0] 

186.08 

280.2679 

Fe^+ 

5.68[1] 

1.119[1] 

9.326[0] 

5.594[0] 

20.14 

257.7285 

Go^°+ 

7.30[1] 

1 . 201 [ 1 ] 

1 . 001 [ 1 ] 

6.006[0] 

-186.87 

238.9746 

Ni“+ 

9.20[1] 

1.261[1] 

1.051[1] 

6.305[0] 

-505.53 

222.5113 

Cul2+ 

1.14 [2] 

1.309[1] 

1.090[1] 

6.543 [0] 

-932.75 

208.1047 

Zni3+ 

1.40 [2] 

1.347[1] 

1.123[1] 

6.736[0] 

-1482.74 

195.3790 


NOTES. -List of the M2 transition rates, Am 2 (in s“^), the reduced transition rate (in 10^^ 

A® • cni~^ • T“^ • s“^), and the wavelength A (in A) of the MIT-transition 3p‘^3d “* 07/2 —t 3p® ^P 3/2 together with 
the energy splitting AE = i5(^D5/2) — E(^D 7 / 2 ) (in cni~^) for Cl-like Ar to Zn. x[n] indicates x x 10'“. 




element (Z) 


Figure 3; (Color online) The total transition rate A = Am 2 + Am it of the 3p^3d '^ 07/2 —>■ ^P 3/2 

transition along the Cl-like isoelectronic sequence for some magnetic-field strengths. We used the fine 
structure energy of (a) 20.14 cm~^ and (b) 5 cm~^ for iron, respectively. 

























7 


Table 3: CONVERGENCE STUDY OF THE CALCULATIONS 


layer 

W 

S 

DF 

0.5305 

4.063 [4] 

n=4 

0.5285 

3.555 [4] 

n=5 

0.5284 

3.355 [4] 

n =6 

0.5283 

3.298 [4] 

n=7 

0.5283 

3.264[4] 


NOTES. -Convergence study of calculated off-diagonal reduced matrix elements W = (^D 5 / 2 ||N*^^^-|-AN*^^^ || '‘D 7 / 2 ) 

I 12 

(in a.u.) of the magnetic interaction and line strength S = (^^ 3 / 2 ||"^D 5 / 2 ) (in a.u.). DF represents a single 
configuration Dirac-Hartree-Fock model. x[n] indicates x x 10". 

B. Fe X 

Due to the close to complete cancellation for iron of the energy difference between the two 
'^D-levels (see Fignre[2]) we will pay special attention to this ion. 

It is clear that some of the properties in Table [2] are more easily obtainable throngh theoretical 
calculations. We illnstrate this in Table [3] where we show the convergence of the calcnlated off- 
diagonal rednced matrix elements, W = (^ 05 ^ 2 11^ 07 / 2 ); representing the magnetic 

interaction, together with the line strength, S = | (^ 7 - 3 / 2 11 1 ^ 05 / 2 ) |^, of the close-lying El tran¬ 
sition. Since these valnes converges fast and are not subjected to cancellation effects, we estimate 
their accuracy to be well within a few percent. This will in turn imply that the prediction of the 
rednced transition rate a^jrp[M) in Equation 0] and [5] is of similar accnracy. 

We give in Eignre 0] the energy difference AE as a fnnction of the maximum n in the active 
set and thereby of the size of the CSE-expansion. It is clear that we also here reach a convergence 
close to a few cm“^ for this property. The final valne for the fine strnctnre splitting is 20.14 cm“^, 
in good agreement with the result from recent confignration interaction calculations (Deb et al. 
2002) as well as Many-Body Pertnrbation Theory (Ishikawa et al. 2010). This strongly snpports 
the prediction of the close degeneracy of the two levels for iron. 

The strong resonance effect for iron is especially fortnnate due to its high abundance in many 
astrophysical plasma. As a matter of fact, the gronnd transition is one of the “coronal” lines used 
to determine the temperatnre of the corona and is known as the corona red line (Swings 1943). 
Unfortnnately there is no firm experimental valne for the critical energy splitting between the 
^D 7/2 and ^D 5 / 2 - At the same time, it is a great challenge to calcnlate the size of this accidental 
degeneracy accurately. 

In the early experimental work by Smitt (Smitt 1977) the two levels were given identical exci¬ 
tation energies of 388708 cm“^. Since then there has been a great deal of work by different groups 
to study the structure of Fe X (see Table 0]) . The differences between the various calculations are 
often mnch larger than the predicted fine-structure splitting, which leads to large nncertainties in 
the level ordering and line identifications. Hnang (Huang et al. 1983) and Dong et al (Dong et al. 
1999) performed mnlti-confignration Dirac-Fock (MCDF) calcnlation and predicted the ^D 5 / 2 , 7/2 
levels to be separated by 78 cm“^ and 108 cm“^ respectively. Ishikawa (Ishikawa et al. 2010) pre¬ 
dicted 18 cm“^ from his Mnltireference-MBPT method. There is a recommended valne from solar 
observations of around 5 cm“^, determined from short-wavelength transitions from higher levels 
and therefore probably qnite nncertain. Predictions from The Goddard Solar Extreme Ultraviolet 
Rocket Telescope and Spectrograph SERTS—89 (Thomas et al. 1994) and SERTS—95 (Brosins et 
al. 1998) spectra give the same energy for ^D 7/2 and ^D 5 / 2 , probably due to limited resolution. 
Einally Del Zanna (Del Zanna et al. 2004) benchmarked the atomic data for Fe X and suggested 
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Figure 4: Convergence trend of the fine-structure separation AE = i?(^D 5 / 2 ) — ^^(^ 07 / 2 ) for Fe®'*' as the 
size of the active set of orbitals is increased as defined by the maximum n-quantum number. 

the best splitting energy to be 5 cm“^. Although our calculations reach a convergence within the 
model, to the final value of around 20 cm“^, it is clear that systematic errors, such as omitted 
contributions to the Hamiltonian, could be relatively important in estimating the accidental de¬ 
generacy of the two levels. We use both our theoretical value and the recommended solar spectral 
value to illustrate the dependence of the average rate Am it on the magnetic field in Figure O It is 
clear that even for relatively weak magnetic fields of only a few hundreds or thousands of Gauss, 
the Amit will be significant compared to the competing M2-rate. 


C. Experimental Determination of The Energy Splitting 

To improve the accuracy of the estimated rate of the MIT, we need to turn to experiment for an 
accurate determination of the energy splitting. For this, we need to overcome two difficulties, first 
enough spectral resolution, and second a light-source with a low electron density and a magnetic 
field. The first requirement is not impossible to fulfill since the fine structure separation can 
be determined using a large spectrometer with a resolution of around 80000. This is far from 
the highest resolution achieved, since e.g. a spectrometer at the Observatory in Meudon has a 
resolution of 150000. However, the line from the ^ 07/2 level has not been observed due to strict 
requirements on the light source. Most sources used at Meudon have generated too dense plasmas 
in which photon transitions from long lived levels cannot be seen (collisions are destroying the 
population of the upper state before the photon is emitted). In addition to this, observation of the 
^D7/2 —P3/2 line requires a strong enough magnetic field of, say a, tenth of a Tesla. This leads 
arguably to only two possible light sources on earth: Tokamaks (Wesson 2004) and the Electron 
Beam Ion Traps(Levine et al. 1988) (EBITs). Tokamak plasma may be too dense, but since 
the magnetic fields involved are higher than what we are discussing here, the line might still be 
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Figure 5: A plot of the Amit as a function of magnetic fields at AA = 20.14 cm ^ and AA = 5 cm and 
compared to Am 2 - 


observable. However, the best choice for our purposes is an EBIT, which has an inherent magnetic 
field to compress the electron beam and is a low density light source. Although the Meudon 
spectrometer demonstrates that the required resolving power can be achieved this instrument is 
not computable with the EBIT operating parameters and a dedicated instrument is required. 

To illustrate the usability of the EBIT source, we have made several model calculations to 
predict the relative strength of the two involved transitions under different circumstances. It 
should be made clear that the EBIT is a light source with a mono-energetic beam of electrons and 
that these models therefore are designed to predict conditions different from those in solar flares 
or the corona. It is important to remember that the intermediate goal before we can proceed is to 
propose an experiment to determine the crucial ^ 07/2 energy separation. We present model 

calculations of the line ratio as a function of magnetic held and electron density (Figure [6]) and 
magnetic held (Figure [7]) of the EBIT, based on collisional-radiative modeling using the Flexible 
Atomic Code (Gu 2008). We show in Figure [6l for several magnetic helds, how the ratio of the 
rates of the magnetic-held induced ^^t /2 P3/2 and the allowed ^ 05/2 P3/2 transitions varies 
as a function of the electron densities. It is clear that the magnetic-held induced line is predicted 
to be visible for the typical range of electron densities of an EBIT, that is 10® — 10^^ cm“® (this 
happens to coincide with the range for solar hares). It is also clear from Figure [71 where we show 
the dependence of this ratio on the magnitude of the external magnetic held for a hxed density of 
10^^ cm~®, that the line ratio will be sensitive to the magnetic held strength. 
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Electron density (cm'^) 

Figure 6 : The ratio of the rates for the magnetically induced ^ 07/2 ^ ^P 3/2 line and the allowed 

^ 05/2 ^ ^P 3/2 line as a function of electron density and magnetic field in an EBIT. Calculations is 

for a monoenergetic beam energy of 250 eV and these data are displayed for some selected magnetic field 
strengths and over a range of densities (10® — 10^^ cm“®) which covers the range for solar flares in the 
corona. Here we used the fine structure energy of 5 cm“^. 

V. CONCLUSION 

To conclude, in this paper we propose a novel and efficient tool to determine magnetic field 
strengths in solar flares. The method is useful for cases of low densities and small external magnetic 
fields (hundreds and thousands of Gauss) that have so far eluded determination. We illustrate 
that a spectral feature originating from the Fe®^ ion is of special interest since it shows a strong 
dependence on the magnetic field strength, with two spectral lines drastically changing their relative 
intensities. We propose a laboratory measurement of the fine structure energy separation between 
the two involved excited states, a crucial parameter in the determination of the external field. When 
this energy separation has been established one can use our theoretical values for the reduced rate of 
the magnetic-field induced transition, which have an accuracy to within a few percent, to calculate 
the atomic response to the external magnetic field. Armed with this it is possible to design a space- 
based mission with a probe that could continuously observe and determine the reclusive magnetic 
fields of the solar flares. 
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Magnetic field (Gauss) 
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